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Q ! Abstract 

^ I The diffuse meta-galactic radiation field at ultraviolet to infrared wavelengths - 

'^ ■ commonly labeled extragalactic background light (EBL) - contains the integrated 

I . emission history of the universe. Difficult to access via direct observations, in- 

^ '. direct constraints on its density can be derived through observations of very-high 

^ , energy (VHE; E>100 GeV) y-rays from distant sources: the VHE photons are at- 

tenuated via pair-production with the low energy photons from the EBL, leaving 
a distinct imprint in the VHE spectra measured on earth. Discoveries made with 
^ ■ current generation VHE observatories like H.E.S.S. and MAGIC enabled strong 

)0 ! constraints on the density of the EBL, especially in the near-infrared. In this ar- 

ticle the prospect of future VHE observatories to derive new constraints on the 
EBL density are discussed. To this end, results from current generation instru- 
VT) • ments will be extrapolated to the future experiment's sensitivity and investigated 

for their power to enable new methods and improved constraints on the EBL den- 
sity. 
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1. Introduction 

The observation of very -high energy y-rays (VHE; E > 100 GeV) from distant 
sources offers the unique possibility to probe the density of the meta-galactic ra- 
diation field at ultraviolet (UV) to infrared (IR) wavelengths, which is commonly 
labeled the extragalactic background light (EBL; typically 0. 1-100/im). The VHE 
y-rays interact with the low energy EBL photons via the pair production process 
(TvheTebl -^ e^e~) and the flux is attenuated III, 120. This attenuation can leave 
distinct signatures in the measured VHE spectra. With assumptions about the 
source physics and the spectrum emitted at the source location (intrinsic spec- 
trum), constraints on the density of the EBL can be derived [e.g. 3, 4]. 

The current generation of VHE instruments (H.E.S.S., MAGIC, VERITAS) 
significantly increased the number of known extragalactic VHE sources from 4 
in the year 2003 to more than 25 today. These discoveries, combined with the 
advanced spectral resolution of these instruments and the wide energy range they 
cover, led to new strong constraints on the EBL density, in particular at optical to 
near-IR (NIR) wavelengths [[Sl, [g, lZD. Since these limits depend on assumptions 
about the source physics, the strong constraints also sparked intense discussions 
on the validity of the assumptions and possible caveats re.g.l8l.l9llld. lllUl2l ll3n. 



This discussion has not yet converged and there are interesting arguments for both 
sides. 

Current generation systems have recently been upgraded (MAGIC-II) or the 
upgrades are under construction (H.E.S.S. II). These upgrades are mainly aimed 
to improve the overall sensitivity by a factor two to three and extend the energy 
range toward the lower energy regime of 20 to 100 GeV. This will lead to some 
improvements, but a quantitative difference or a breakthrough compared to the 
performance of the existing facilities will only be achieved with an order of mag- 
nitude improvement in sensitivity. The Next Generation Cherenkov Telescope 
Systems (NGCTS) are in the advanced planing phase aiming to achieve this order 
of magnitude improvement: the Cherenkov Telescope Array (CTAIj) lll4ll and the 
Advanced Gamma-ray Imaging System (AGIScI) lllSn . Whereas CTA envisions 
to improve the sensitivity over a wide energy range from the few tens of GeV to 
the multi TeV regime, AGIS mainly concentrates on energies above 100 GeV, an 
extended field of view and an improvement of the angular resolution. 

The potential of these upcoming experiments to probe the EBL is the topic of 



'http://www.cta-observatory.org/ 
^http://www.agis-observatory.org/ 
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Figure 1 : Present day (z - 0) EBL density of the EBL models/shapes utilized in this study (Stecker 
et. al: [16]: Raue & Mazin: 1171 : Franceschini et al. Iil8i1 ). Grey markers show measurements and 
limits on the EBL density (from ||6[]). 

this article. While an order of magnitude improvement in sensitivity for an as- 
tronomical instrument will always lead to new and unexpected results, this article 
will - as a first step - focus on known results and their extrapolation according 
to the sensitivities of the next generation instruments. Emphasis will be on new 
techniques enabled by the performance features (extended sensitivity and energy 
range) of the upcoming instruments. 

For the calculations in the paper a standard ACDM cosmology with h = Qa = 
0.7 and Q.m = 0.3 is adopted. 



2. Basic assumptions and simulation details 

EBL models and attenuation. The precise level of the EBL density is not well 
known. Solid lower limits from integrated deep galaxies counts at optical and 
infrared wavelengths do exist [e.g. Il9l. l20l l2lL 12211 . but direct measurements are 
hampered by dominant foregrounds [23] . Over a large wavelength region in the 
infrared the best upper limits on the EBL density are derived from VHE observa- 
tions of distant sources [e.g. a [g, lZ|]. To account for the uncertainty in the EBL 
density two different approaches are followed here: (i) To illustrate the effect of 
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Figure 2: Attenuation for VHE y-rays for the EBL models utilized in this study for sources at 
redshift z - 0.03 (thin lines) and z - 0.2 (thick hnes). 



the different EBL densities two extreme EBL models are used: for the low EBL 
density case the model from [18] (FranOS in the following) is adopted, for the high 
model the fast evolution case from [16] (SteckerFE in the following) is used. It 
should be noted that the later model is already disfavored by several VHE obser- 
vations, (ii) Detailed studies of the effect of different EBL densities are carried out 
by scaling the EBL density presented in [17] (Gen in the following) and FranOS. 
The present day EBL densities (z = 0) for the EBL models utilized in this work 
are displayed in Fig. \T\ 

Fig. [2] displays the resulting attenuation for VHE y-ray sources at 2 differ- 
ent redshifts (z = 0.03 and 0.2) for the models utilized. Several features can be 
identified: 

• At low energies (< SOGeV) the spectrum is practically not attenuated. 
Since this energy range will be sampled by the NGCTSs with high pre- 
cision, it will be possible to measure the unabsorbed spectrum. 

• At energies between SOGeV and 2 TeV the attenuation is increasing due 
to the EBL photons in the optical to near-infrared range peak of the EBL 
density. 



• At energies between 2 to 10 TeV a flattening of the attenuation is expected, 
due to the ~ A'^ behavior of the EBL density in the near to mid-infrared, 
resulting in a constant attenuation. Such a modulation of the EBL attenu- 
ation has been considered as a possible key signature for EBL attenuation 



[e.g.|24j,|25|]. Unfortunately, the intrinsic weakness of the sources combined 
with the sensitivity of the instruments make it very difficult to probe such 
a feature with previous or current generation experiments. For the Steck- 
erFE model the EBL density in this wavelength range is flatter, resulting in 
a smoother attenuation from 100 GeV to 10 TeV, suppressing such a feature. 

• At energies around 10 TeV the turnover in the EBL density towards the far- 
infrared peak of the EBL results in a strong attenuation, effectively resulting 
in a cut-off in the measured spectra. 

The strength and the position of these features vary with the distance of the VHE 
sources, the assumed EBL model, and the overall EBL density. 

EBL limits from VHE observations. So far, VHE sources used to derive limits on 
the EBL density belong to a single source class, active galactic nuclei (AGNs), 
and the majority of them to the Blazar sub-class, which are AGNs with strong jet 
activity and the jets are closely aligned to the line of sight of the observer. Up 
to now, mainly two different methods - and thereby assumptions about the source 
intrinsic spectrum - have been utilized to derive limits on the EBL density: 

• Spectral concavity. It is assumed that the overall intrinsic source spectrum 
at high energies will follow a concave shape, or at least will not show an 
exponential rise towards the highest energies. These assumptions are well 
motivated by the common leptonic modeling of the sources under investiga- 
tion (blazars), although different (maybe more exotic models) can possibly 
reproduce such a feature [e.g.|8|]. Limits on the EBL density are derived by 
excluding EBL densities that would lead to such features in the observed 
sources. This method naturally probes the EBL at wavelengths from the 
mid to the far-infrared. 

• Maximum spectral hardness. To probe the EBL in the optical to near- 
infrared, it is assumed that the intrinsic source spectrum cannot exceed 
a certain absolute hardness. While somewhat similar in spirit to the first 
method the underlying assumptions are stronger, since in the energy range 
of interest (100 GeV to several TeV) the spectral shape of the intrinsic spec- 
trum is more uncertain. While most of the basic models used to describe the 
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Figure 3: Integral flux sensitivity (5cr in 50 h) of the NGCTS used in this study in com- 
parison to the sensitivity of existing observatories (H.E.S.S.: ll26ll : MAGIC: 112711 : Fermi: 
http : //www-glast . slac . Stanford. edu/so£tware/IS/glastJ.at_Performance . htm). 



source spectra indeed imply that the VHE spectrum does not exceed a cer- 
tain hardness, the absolute value is less certairO and possible source intrinsic 
effect (e.g. internal absorption [11]) could complicated the situation. 

In this paper two different methods to derive limits on the EBL density will be 
explored: (i) utilizing the unabsorbed part of the VHE spectrum and (ii) searching 
for attenuation modulation signatures. While not completely new, it will be shown 
that with the NGCTS 's extended energy range paired with its vastly improved 
sensitivity it will be possible to utilize these methods effectively for the first time. 
Method (i) holds the potential to derive limits on the EBL density with a minimal 
set of assumptions, while method (ii) enables to not only derive upper limits on 
the EBL density but to probe the absolute level. 

Effective area, background rate and sensitivity. To simulate the spectra, which 
will be detected by the NGCTS, assumptions about the sensitivity of the instru- 



e.g. r = 1.5 - 0.6 for dN/dE ~ E ^ and simple leptonic models. 



ment have to be made. As baseline sensitivity, the sensitivity of the "4 large + 



85" CTA array presented in [12611 is taken and the sensitivity is derived assuming 
an effective detector area (effective area) and a background rate. For the effective 
area the post cut MAGIC effective area at 20deg zenith angle ( [.28.1 ) scaled up 
by a factor 20 is adopted to reach an effective area exceeding lO^m^ at energies 
above a few hundred GeV. In addition, the MAGIC effective area is shifted by a 
factor of 2 to lower energies to reflect the improved sensitivity at low energies. □ 
The resulting effective area Angcts versus energy E in the energy range 20 GeV 
to 20 TeV is well described by the function 

Angcts = 2.5 ■ 10^W{E/l TeV)^(l + (£/(l TeV ■ Om)f'^r^'°'^ (1) 

The differential background rate after event selection cuts is approximated by a 
broken power law function with photon index 3.6 below and 2.7 above the break 
energy of 300 GeV, similar to what is seen in current generation instruments. The 
absolute level of the background flux rate is chosen so that the resulting sensitivity 



matches the sensitivity of the CTA "4 large + 85" array from [12611 . The integral 
sensitivity adopted in this study for the NGCTS in comparison to current genera- 
tion and future HE/VHE instruments is displayed in Fig. |3l In the overlap re gion 
it follows very well the CTA sensitivity for the "4 large + 85" array from [261 
up to 10 TeV. Compared to current instruments (H.E.S.S., MAGIC, VERITAS) 
the sensitivity in the core energy range between 100 GeV and 10 TeV is improved 
by about one order of magnitude. In addition, the energy range is extended to- 
ward lower and higher energies. With the chosen parameters for the effective area 
and the background rate at low energies significant sensitivity is reached down 
to energies below 20 GeV, which enables a large overlap region in energy with 
the Large Area Telescope (LAT) onboard Fermi satellite. For this energy region, 
where there is no overlap with the published CTA sensitivity (<40GeV), there is 
of course a certain degree of freedom in the choice of parameters and therefore 
the sensitivity is not very well constrained. In this study only energies between 
40 GeV and 10 TeV wiU be used. 

Spectrum simulation method. To calculate the simulated spectrum the number of 
y-photon events Ns (signal) and background events Nbg per energy bin need to be 



"*The energy threshold approximately scales linearly with 1/Area of the mirror area of the single 
telescope. For a CTA type NGCTS instruments the largest telescopes are expected to have mirror 
diameter in the order of 24 m which gives a factor two larger area compared to 1 7m for the MAGIC 
mirror 



determined: the NGCTS effective area folded with the assumed intrinsic source 
flux is integrated over the energy bin; the same is done for the background events 
integrating over the background rate in the bin. Both numbers are multiplied 
with the eff"ective observation time. Diff"erent functions are utilized to describe 
the intrinsic flux and they will be discussed further in the section where they first 
appear. The attenuation of the source flux due to the EBL is calculated following 
the recipe given in [|6L HSfl: the attenuation is directly folded into the intrinsic 
spectrum and then the attenuated intrinsic flux function folded with the efi'ective 
area is integrated over the energy bin. The number of signal events is randomly 
varied assuming a Poisson distribution. It is assumed that the background is well 
determined (e.g. via background measurements in a large sky area compared to 
the signal region). The error on the number of photons A^^- in a bin is derived 
utilizing equation 17 from [i29ll assuming 5 background regions (i.e. an alpha 
factor of 0.2). The signal in an energy bin is considered significant when all 
of the following criteria are met: (1) the signal significance exceeds 3 standard 
deviations, (2) there are at least 10 excess events in the bin, (3) the number of 
excess events in the bin exceeds 3% of the number of background events. This last 
condition takes into account a systematic error in the determination of the number 
of background events. These are rather conservative assumptions, since e.g. in 
current publications on VHE y-astronomy often energy bins with significances 
down to 1.5 cr or less are included in the analysis. 

lACT experiments have a limited energy resolution which, for current gener- 



ation instruments, is in the order of <15% for energies above 100 GeV [e.g. |30D 
and <40% down to 70 GeV [28]. For MAGIC-II an energy resolution of <25% 
for energies down to ~ 50 GeV is achievable [see Fig. 3 of 31]. Due to the, on 
average, higher number of telescopes participating in each event and the increased 
mirror size of the large telescopes the energy resolution for an NGCTS is expected 
to improve further. Such an energy resolution paired with an energy spectrum un- 
folding method [e.g. '321] will enable to robustly reconstruct smooth spectral shapes 
(e.g. power law or log parabola) even down to low energies. The reconstruction 
of an EBL attenuation structure at mid-energies is possibly more aff'ected by the 
limited energy resolution and this will be further discussed at the end of Sect. 4. 

Simulation example. Fig. |4] displays simulated spectra for two sources and two 
assumed EBL densities for an observation time of 20 h. The intrinsic spectrum 
is assumed to follow a simple power law (0(£') = Oq x {Ejl TeV)"'^), with the 
parameters adopted so that the EBL attenuated simulated spectrum matches the 
measured one. Shown are results for lES 1101-232 {z = 0.186), a hard spectrum 
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Figure 4: Simulated VHE spectra for two sources and two EBL models (blue:Fran08; 
red:SteckerFE). Upper panel: IBS 1101-232. Lower panel: PKS2155-304. Grey markers show 
the measured spectra and the measured spectra de-absorbed for EBL attenuation. Red and blue 
lines are the assumed intrinsic spectra emitted at the source for two different levels of the EBL 
density, matching the measured de-absorbed spectra. Red and blue markers show the expected 
spectrum as measured for a next generation VHE instruments simulated in this work. The dashed 
line gives the flux of the Crab Nebula. Spectral points below 40 Ge V are only shown for illustrative 
purpose and are not used in the analysis. 



i^vHE ~ 3) distant source, whose discovery at VHE energies enabled to derive 
strong limits on the EBL density in the optical to near-infrared [5], and PKS 2155- 
304, which is a fairly strong VHE source (about 20% Crab in the quiescente state) 
with a softer spectrum {Tvhe ~ 3.3) located at an intermediate redshift of z = 
0.116. In the case of lES 1101-232 the hard measured spectrum in combination 
with the larger distance leads to a very hard intrinsic spectrum (i.e. harder than 
anticipated by simple leptonic and hadronic models), even for the low EBL model. 
Such a hard intrinsic spectrum coupled with the relative weakness of the source 
(~2% Crab) leads to a simulated spectrum which does not largely increase the 
energy range covered compared to the H.E.S.S. measurement of lES 1101-232: 
neither at low energies (due to the hard spectrum) nor at high energies (the strong 
EBL attenuation suppresses the signal below the NGCTS sensitivity). Even for 
the case of the two very different assumed EBL models, no strong difference in 
the simulated spectra is apparent, which could be used to differentiate between 
the models. This is different in the case of PKS 2155-304, where, for the different 
EBL models, very different spectra are expected to be measured at lower energies. 
The behavior at low energies will be discussed further in the next section. 

3. Utilizing the unabsorbed part of the spectrum 

One of the main features of a NGCTS will be a high sensitivity in the energy 
range between 20 and 100 GeV, an energy range which holds the possibility to 
directly sample parts of the energy spectrum of a source, which are not affected 
by the EBL attenuation. In this section it will be explored how this energy range 
can be utilized to derive limits on the EBL density, and what are possible problems 
and caveats. 

Simulation & analysis chain. The following simulation and analysis chain is uti- 
lized: 

1 . Calculate EBL attenuation for a specific EBL density and source distance. 

2. De-attenuate a measured spectrum and fit with power law {dN/dE = Oq • 
E~^). The fit results serve as input flux function for the simulated spectrum. 

3. Simulate spectrum as measured by an NGCTS with calculated input flux 
function. 

4. Fit simulated spectrum in a low energy regime (intrinsic spectrum) and high 
energy regime (absorbed spectrum). Again, a simple power law function is 
used in each energy regime. 
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Figure 5: Example of simulated spectra for different EBL densities. The base spectrum assumed 
is the quiescent state spectrum of PKS 2155-304 (z = 0.1 16), EBL model is from FranOS. Shown 
are: the measured spectrum (grey markers), the simulated spectra for different level of the EBL 
density (black markers) and the corresponding assumed intrinsic spectra (black lines), the source 
spectrum in the GeV energy range as measured by Fermi |33] (purple butterfly), and the energy 
ranges which are used to determine the slope of the simulated spectrum at low (blue) and high 
(red) energies. Spectral points below 40 GeV are only shown for illustrative purpose and are not 
used in the analysis. 
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Figure 6: EBL density (z - 0) for the FranOS EBL model scaled with a factor 0.7 to 1.7 in steps of 
0. 1 . The red dashed curve is the unsealed EBL model. 



5. 200 spectra are simulated and mean values are used. 

An example for this procedure is shown in Fig. [5] for the quiescence spectrum of 
PKS 2155-304 (z = 0.1 16). An observation time of 20 h has been assumed, which 
could easily be extended given that this is a steady flux state. As EBL model the 
FranOS model is used, scaled in steps of 0. 1 from 0.7 to 1 .7 (Fig. (6]). 

Fit-range &- source intrinsic break. The resulting EBL attenuations for the scaled 
models are shown in Fig. |7] upper left panel. The energy ranges used to fit the 
power laws are marked by colored boxes. In the high energy range the attenu- 
ation follows approximately a power law and since the input spectrum is also a 
power law, the resulting attenuated spectrum will follow a power law as well. In 
the low energy range no significant absorption is present. As recent Fermi/LAT 
observations show, a spectral break is observed between the GeV and the TeV 
energy range [|34n . For many sources this break can be attributed to a break ex- 
pected from EBL attenuation. For these sources an NGCTS would be able to 
sample the intrinsic spectrum down to very low energies. In other cases, e.g. 
PKS 2155-304, the break between GeV and TeV range is stronger than expected 
from EBL attenuation, therefore an additional, source intrinsic break is expected 
somewhere between the two energy bands. Currently, the statistics for this energy 
range for the sources considered here is not yet sufficient to correctly model the 
break. Therefore, for the analysis in this paper the lower edge of the low energy 
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Figure 7: Results from the power law fit to the simulated spectrum in the low and high energy 
band using the PKS 2155-304 quiescence spectrum with the FranOS (upper row) and Gen (lower 
row) EBL model. Left: Attenuation of the VHE y-ray flux resulting from the scaled EBL models. 
The colored boxes in the background mark the fit range for the power law at low (red) and high 
(blue) energies. Middle: Spectral index F from the fit of a power law to the low (red) and high 
(blue) energy range versus the scaling factor of the EBL model. The color shaded bands denote 
the error on the spectral index from the fit (RMS of the spectral index distribution). Black crosses 
mark the intrinsic spectral index that has been utilized. Right: Spectral break between low and 
high energies (error bands from error propagation). For discussion see main text. 
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fit range is chosen to start above the Fermi/LAT energy band (defined as the 5th 
highest photon in energy reported in [33], e.g. ~40 GeV for PKS 2155-304). Fu- 
ture observations with Fermi/LAT, MAGIC-II, and H.E.S.S. II will provide further 
information on this issue. 

EBL attenuation spectral break. Fig. |7] upper middle panel shows the spectral 
index resulting from the power law fit for the two energy bands for the different 
scalings of the EBL model. The markers show the mean spectral index of the fits, 
with the error given as the RMS of the mean spectral index distribution (shaded 
bands). As alternative error definition the mean error on the spectral index from 
the fit could be used which gives similar results. The black crosses mark the 
spectral index utilized for the input source spectrum. It can be seen that the fit in 
the low energy band reproduces very well the assumed intrinsic spectral indices c| 
For the highest scaled EBL densities the effect of attenuation becomes relevant 
in the low energy band and the spectral index from the fit is steeper (larger) than 
the one from the input spectrum. This effect will be discussed in more detail 
below. The right panel of the figures displays the strength of the spectral break 
between the two power laws Ylqw - ^high, with the error bands calculated via 
error propagation. An EBL density scale factor of 0. 1 corresponds to one to two 
standard deviations difference in the strength of the break, u Note that a scaling 
of 0.1 corresponds to an EBL density of ~1 nWm"^sr"^ at 2yum which is of the 
same order as the error on the lower limits from integrated source counts at this 
wavelength. The detection of a break with certain strength can be converted into 
an upper limit of the EBL density (thick dashed line). In principal, a single well 
measured spectrum is sufficient to derive such limits. Since EBL attenuation is a 
global, redshift dependent phenomenon, in general, a combined fit to all available 
VHE data should be used to derive limits on the EBL density [see e.g. u\. This is 
discussed in more details in Sect. [51 

EBL attenuation at low energies. The same analysis has been performed utiliz- 
ing the Gen EBL model and the results are displayed in Fig. |7] lower row. The 
Gen EBL density in the ultraviolet to optical is higher than in the FranOS model 
(Fig.[T]), therefore resulting in a non negligible attenuation in the low energy band 



""The correlation between the EBL density and the spectral index in the low energy band is a 
result of the methodology i.e. for each EBL density a corresponding intrinsic spectrum is con- 
structed. 

^I.e. if the FranOS EBL model is correct and the intrinsic spectrum follows a power law in the 
energy range considered a L3 scaling of the model could be excluded with 3 standard deviations. 
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Figure 8 : Energy dependent upper limits on the attenuation for the FranOS (left) and the Gen (right) 
EBL models derived utilizing the PKS 2155-304 quiescence spectrum. Shown are the limits taking 
into account the 1 cr error on the flux points (orange filled squares), with the addition of the 1 cr error 
on the spectral index of the low energy power law (red filled triangles), and no errors (purple filled 
circles). The black line shows the attenuation for the EBL model utilized, grey lines show the 
attenuation for the EBL model scaled in steps of 0. L 



(Fig. H lower left panel). Consequently, the power law fit in the low energy band 
results in a too soft spectrum compared to the input source spectrum (Fig.|7]lower 
middle panel). This mis-reconstruction, if taken at face value, would lead to an 
overestimated upper limit on the EBL density (Fig. |7] lower right panel). To be 
able to utilize the low energy part of the spectrum as proxy for the intrinsic spec- 
trum it is therefore crucial to carefully examine the spectral shape in the energy 
range for signs of curvature. A conservative approach would be to correct the flux 
points in the fit energy range for the attenuation from a maximum EBL before 
performing the fit. On the other hand, if a high statistic VHE spectrum from a 
distant source does not show any indication of curvature in this energy regime, 
tight limits on the EBL density in the ultraviolet to optical EBL can be derived. 

Energy dependent attenuation limits. The power law fit to the low energy region, 
if taken as proxy for the intrinsic spectrum, can also be used to derive energy 
dependent upper limits on the attenuation. Such limits can be extremely useful: 
the EBL attenuation is an energy dependent process, where certain wavelength 
regions of the EBL spectrum are connected to certain VHE energy ranges. Fea- 
tures in the EBL density, as e.g. expected from the first stars [see e.g. |35|], could 



15 



therefore also produce features in the attenuation. The resulting energy dependent 
EBL attenuation limits from such an analysis for the FranOS and the Gen EBL 
model utilizing the PKS 2155-304 quiesent spectrum are shown in Fig.HQ. As can 
be seen from the figure, the dominant source of errors in such an analysis are the 
uncertainties on the spectral index from the power law fit in the low energy range, 
since the statistical errors on individual flux points are comparably small. Again, 
the rejection power is in the order of 1-2 cr per 0.1 scaling of the EBL density at 
energies > 1 TeV, when utilizing a simple and ad-hoc error definition. In the case 
of the Gen EBL model again the effect of the attenuation "spill over" in the low 
energy fitting band is visible: due to the EBL attenuation the power law fit is too 
soft, resulting in an underestimated EBL attenuation and too strong limits. 

Flaring states. Up to here only the quiescent flux state of sources has been con- 
sidered, since its observation is guaranteed. AGNs do also show extreme flaring 
behavior with high fluxes, unfortunately usually only for short periods of time. 
Here, the VHE spectrum recorded during the extreme flux outburst of PKS 2155- 
304 in 2006 [36] is investigated. The spectrum is an averaged spectrum of 1 .5 h 
of observations, in which the source showed several strong flares with time-scales 
down to a few minutes. While no strong spectral variations in the VHE spectrum 
have been observed during this flare, for the analysis only short observations of 
6 min length have been considered, demonstrating the power of the NGCTS to 
derive high quality spectra even on such short time- scales. In addition, the EBL 
attenuation is a "stationary" effect in the spectrum, therefore a temporally fine 
resolved flaring state might enable to distinguish between intrinsic and EBL ef- 
fects. The resolution of the EBL density achieved for such a 6 min observation 
of the source in a flaring state is comparable with the one achieved for 20 h of 
quiescence observations (Fig.|9]top row). 

Nearby sources. Strong flux outbursts in the VHE regime have also been observed 
in the two nearby sources Mkn421 and Mkn501. While the smaller distances 
result in less EBL attenuation and therefore a weaker signature, the overall higher 
flux levels of these two sources enable to derive spectra with high event statistics. 
Here, the time-averaged VHE spectrum of Mkn501 during the 1997 high state is 
investigated, assuming an observation time of 20 h. As can be seen from Fig. |9] 
bottom row, while the overall strength of the break between the two power laws is 



^For this investigation only a single simulated spectrum without applying the statistical dicing 
of the event numbers is used. 
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Figure 9: EBL limit results from power law fits to low and high energy bands of the VHE spectrum 
using scaled version of the FranOS EBL model (see Fig. [T] and [8] for detailed description). Upper 
row: Results for the VHE spectrum from the flaring state of PKS 2155-304 in 2006 |36] (6min 
observation time), Lower row: Results for the averaged high state spectrum of Mkn501 1371] 
(z = 0.034; 20 h observation time). 
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much smaller than in the case of PKS 2155-304 (0.15-0.4 vs 0.5-1.2) comparable 
limits can be derived due to the superior statistics. In addition, the problem of 
attenuation in the low energy band is less severe. One caveat - at least for the case 
of Mkn421 - is the fact that there are strong indications for an intrinsic (not EBL 
related) break between the GeV and TeV range at relatively high energies (around 
100 GeV) [I34I1 . which would make it difficult to define a proper low energy region 
void of any curvature. A different method, discussing the possibilities of how such 
high quality spectra with a wide energy range produced by a nearby source can be 
used to test the EBL density, is discussed in the next section. 

4. Attenuation modulation at mid-energies 

The NGCTS will achieve an unprecedented sensitivity in the intermediate en- 
ergy range, i.e., between 100 GeV and few TeV. For nearby sources, this energy 
range is most sensitive for the EBL density in the optical to infrared regime (~1- 
15 //m), which can, thus, be probed very efficiently with an NGCTS measurement 
of an AGN energy spectrum. 

Smoothness of AGN spectra. The measured energy spectra of AGNs in the energy 
range between 100 GeV and few TeV follow usually a smooth shape. For most 
of the measured sources, a simple power law fit is sufficient to describe the avail- 
able data well, whereas for sources in a flare state (like the flare of PKS 2155-304 
in 2006) or with a generally high emission state (like Mkn421), either a curved 
power law or a power law with a cut-off are successfully used. The curved power 
law (also known in the literature as the double-log parabola) is expected to de- 
scribe the spectra well at energies close to the position of the Inverse-Compton 
peak. The power law with a cut-off instead is the expected behavior of a source 
which does not provide necessary conditions for acceleration of charged parti- 
cles to sufficiently high energies. All scenarios do have one common feature: the 
measured spectra can be described by smooth functions,i.e., no features, wiggles 
or pile-ups are expected, especially after de-convolving the spectra for the effect 
of the EBL absorption. This property can, therefore, be used to distinguish be- 
tween different overall EBL levels in the optical to infrared regime: whereas the 
"correct" EBL model and level will produce a smooth intrinsic AGN spectrum, an 
"incorrect" EBL level would result in a signature (in form of well defined wiggles) 
in the reconstructed intrinsic spectrum. 

The strength of the method is that the EBL signatures in the reconstructed 
AGN spectra will not only be visible (measurable) in the case where the assumed 
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EBL level is higher than the real one, but also in case the assumed EBL level is 
lower than the real one. It is, therefore, the first indirect method to really measure 
the EBL density at z=0. 

Simulation & analysis chain. The utilized methodology is sketched below: 

1. Assume an intrinsic spectral shape and the flux level of a known strong 
extragalactic gamma-ray source. 

2. Simulate NGCTS spectrum assuming the absorption due to the standard 
EBL, i.e., with the scaling factor of 1. 

3. Reconstruct intrinsic spectrum of the source assuming a scaled EBL level. 
For a scaling, which is diff"erent enough from the standard EBL, the recon- 
structed AGN spectrum will show distinct wiggles. 

4. To characterize the presence of the wiggles, a fit by a smooth source func- 
tion (spectral shape) is performed. The chosen fit shape is the curved power 
law: dN/dE = NoX [E'"^^ iog(£/£o)J 

5. The resulting x^ of the fit is then used to judge if the change in the EBL 
level results in an improbable reconstructed intrinsic spectrum. 

6. Repeat the simulation 1000 times for the given EBL scaling and compute 
the mean and the RMS of ^^ values from the fits to the reconstructed intrin- 
sic spectrum. 

Simulation example. The steps 1-4 of the method are illustrated in Fig. [10]for the 
VHE spectrum of Mkn 501 . The assumed spectral shape and the flux level of the 
intrinsic spectrum are adapted to the flux measured by HEGRA [37] during the 
outburst of the source in 1997: the original data are shown by grey open squares 
in the upper panel of the figure. The simulated NGCTS spectrum calculated using 
the 'correct EBL' (Gen EBL model; scale factor 1) is shown in red, whereas the 
assumed intrinsic spectrum of Mkn 501 is shown by the solid grey line and the re- 
constructed intrinsic spectrum is shown by the blue filled circles. The effect of the 
mis-reconstruction of the intrinsic spectrum is shown for the example of an EBL 
scaled by a factor of 1.3: the reconstructed intrinsic spectrum (green filled trian- 
gles) clearly shows wiggles in the fit range. The effect of the wiggles is more vis- 
ible in the lower panel of the figure where the residuals to the best fit function are 
shown. The wiggles are quantified by a fit in the energy range between 100 GeV 
and 7TeV, well before a possible pile- up in the spectrum arises. The choice of 
the fit range is made in order not to bias the result by the level of the EBL above 
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Figure 10: Search for signatures at mid-VHE in the Mkn 501 spectrum. Top panel: The Mkn 501 
spectral energy distribution (SED) at very high energies. Shown are: measured spectrum by 
HEGRA (grey open squares), simulated NGCTS spectrum (red filled squares), assumed intrin- 
sic spectrum of Mkn 501 (grey solid line), reconstructed intrinsic Mkn 501 spectrum (blue filled 
circles) and reconstructed intrinsic spectrum of Mkn 501 assuming an EBL scaling factor of 1.3 
(green triangles). For comparison, the deabsorbed HEGRA Mkn 501 is shown for the same two 
EBL models (filled grey circles and open grey circles, respectively). The Crab Nebula spectrum 
(grey dashed line) is shown for comparison. Middle upper panel: Residuals of the fits to the 
de-absorbed HEGRA Mrk501 spectrum for the case of the two different assumptions about the 
EBL density. The difference between the two curves (filled and open circles) is very small. The 
HEGRA sensitivity is not sufficient to decide between the two assumptions. Middle lower panel: 
Residuals between the best fit to the SED in case of NGCTS measurement and the spectral points 
from the intrinsic spectrum using the correct EBL density. Bottom panel: Residuals between the 
best fit to the SED for a NGCTS measurement and the spectral points from the intrinsic spec- 
trum reconstructed using the scaled EBL model. A clear and significant signature (wiggles) in the 
residuals is visible, which is quantified by a low probability of the fit. 
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Figure 1 1 : Quantitative results from the search for EBL signatures in the mid-VHE using the 
energy spectrum of MknSOl. Shown are the reduced ;if^ values from the fits to the reconstructed 
spectra of Mkn 501 as a function of EBL scaling factor. Blue filled squares and green filled circles 
show the expected result for 20 and 50 hours of NGTCS observations, respectively. The black 
horizontal lines correspond to fit probabilities as labeled. 



lOfxm, to which the VHE spectra are very sensitive due to a super exponential de- 
pendency of the attenuation with the wavelength in that range. Using the correct 
EBL level to reconstruct the intrinsic spectrum, the intrinsic spectrum is well de- 
scribed by a smooth function (Fig.fTOl middle lower panel). Instead, when using 
a "wrong" scaled EBL density characteristic deviations (wiggles) from a smooth 
function are visible (Fig.[10l lower panel). The wiggles are quantified by a reduced 
X^ value of 51.4/12, corresponding to a fit probability of 7.9 • 10"^. The small fit 
probability for the assumed scaled EBL level implies (under used assumptions) 
significant presence of the unphysical wiggles and, therefore, an exclusion of that 
particular EBL realization. For comparison, HEGRA measured spectrum is also 
de-absorbed using the same two EBL models and the results are shown by grey 
filled and open circles in the upper panel. The residuals to the best fits are shown 
in the middle upper panel. The difference between the two curves (represented by 
grey filled and open circles) is not significant and is difficult to see. The similar 
reduced ;^^ values of 13.6/8 and 14.0/8 underline that the HEGRA measurement 
was not sensitive to the method. 

Results of the analysis. This procedure is repeated 1000 times for every scaled 
EBL density in order to achieve a solid statistical mean and a 1 cr (68%) cover- 
age of the reduced ;^f^ values. The results are shown in Fig. [TT] for two different 
NGCTS exposures (20 and 50 hours) plotting the mean reduced x^ values in- 
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eluding 68% error bars versus the EBL scaling factor. As a result one can see a 
parabola- like curve with the minimum at EBL scaling factor of 1 .0, which is via 
construction the correct EBL model. The expected mean reduced x^ values are 
shown by the blue filled squares and the green filled circles for the exposures of 
20 hr and 50 hr, respectively. The corresponding fit probabilities of P=50%, 1% 
and 0.01% are shown by the dotted, short-dashed and long-dashed lines, respec- 
tively. An EBL scaling is considered to be excluded when the mean reduced x^ 
value including its 68% error exceeds the x^ value for P=0.01%. For the case 
constructed here this means that in case of 50 hours observation, the EBL scalings 
below 0.65 and above 1.35 are excluded. In the case of 20 hours observation, the 
EBL scalings above 1.50 are excluded. 

Caveats of the method. The method presented above is not only sensitive in con- 
straining the EBL density but it is also a first attempt to resolve the actual EBL 
level. Still, the method has several caveats: (i) There is a possibility that the in- 
trinsic VHE spectrum is not smooth for some AGNs. For example, several emis- 
sion regions or several generations of charged emitting particles (e.g., electrons) 
may produce spectral features which would mimic the signature from a wrong 
EBL model. Though neither the AGNs observed at low energies with Fermi/LAT, 
nor the extragalactic VHE y-ray sources observed so far with lACTs have shown 
structured spectra, such wiggles can be interpreted as the result of a wrong EBL 
model only if found to exist in several spectra of different sources consistently, 
(ii) The spectral signature shown in Fig. [10] can only be recognized well after a 
long exposure of a source being in a flare state. The data of the Mkn501 flare 
from 1997 used for the study remains unique so far. This means that it cannot 
be expected that a study with such precision can be done on many sources with 
future NGCTS. Still, the indirect EBL measurement based on few strong flares 
will be of high importance, (iii) The limited energy resolution of lACTs (~10%) 
could affect the possibility of detecting such a signature. The overall shape of 
the attenuation effect is broad enough in energy to not strongly be affected by the 
smoothing of an energy resolution of < 10%. For nearby sources, however, the 
wiggle structure produced by the wrongly assumed EBL model is only of order of 
a few per cent of the total flux level. An additional spread introduced by the spill- 
over effects due to events with mis-reconstructed energy may enhance or alter the 
wiggles. Therefore, an efficient and precise spectral unfolding technique will be 
key to overcome these problems for nearby sources. For more distant sources the 
attenuation signatures get quickly stronger. It will be, therefore, highly interesting 
to study a strong flaring state from a source at redshift > 0.1 with the NGCTS. 
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5. Summary & Conclusions 

In this paper the potential of a Next Generation Cherenkov Telescope Sys- 
tem (NGCTS) to study the EBL through observations of VHE spectra from dis- 
tant sources is explored. In the focus of the study lies the energy range be- 
tween 40GeV and lOTeV, where a factor 10 improvement in sensitivity over 
current generation experiments is expected. Two different methods are investi- 
gated: (i) utilizing the unabsorbed part of the VHE spectrum in the energy range 
40-lOOGeV and (ii) searching for attenuation modulation signatures at energies 
between 100 GeV to 7TeV. While some caveats, like e.g. the exact shape of the 
intrinsic spectra, do exist, overall the two methods show promising results, clearly 
go beyond what is possible with current generation instruments. 

EBL attenuation is a global phenomenon that affects the spectra of all sources 
in the same way. It is therefore quite natural to expand the studies on the EBL by 
combining the results from different methods, sources or source flux states. This 
can e.g. be done by utilizing log-likelihood methods to combine the constraints 
from different sources and methods. In addition only the low and mid-VHEs have 
been investigated. At energies > lOTeV further signatures from EBL attenuation 
("cut-ofF') are expected, which have been used extensively in the past to derive 
limits on the EBL density in the FIR. 

In this paper only data from a single instrument - a NGCTS - and minimal 
assumptions about intrinsic spectrum have been used to derive constraints on the 
EBL density. AGNs are observed all across the electromagnetic spectrum and 
the theory predicts connections between the observations at different energies. A 
wealth of multi-wavelength information is therefore accessible, which - in com- 
bination with a theoretical model - can be used to constrain the spectrum at VHE. 
A first model-dependent approach to determine the EBL level using a consistent 



SED modeling of detected blazars was discussed recently by 113 811 . While a de- 
tailed discussion of this topic is beyond the scope of this paper, as an example it 
should only be mentioned the power of combining the Fermi observations at MeV 
to GeV energies with observations from an NGCTS: such observations will cover 
over 9 decades (!) in energy with high precision. 

The sources under study show strong flux variability in the VHE band. The 
high sensitivity of a NGCTS will enable to study the energy spectrum on very 
short time-scales in great detail. The time resolution will enable stronger con- 
straints on the theoretical modeling of the source spectra since the experimental 
data will be available on shorter scales than the relevant changes in the emission 
regions. This will further help to disentangle source intrinsic effects and EBL 
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attenuation. 

In this study only the VHE spectra of a few selected known sources have been 
investigated. The NGCTS will bring new discoveries. With the high sensitivity in 
the 20 to 100 GeV range it will be possible to detect fainter sources and sources 
at higher redshift. A (sufficiently) large sample of sources will serve a two-fold 
purpose: (1) help to improve the understanding of the intrinsic source physics and 
(2) enable statistical studies of the EBL effects, e.g. by investigating the EBL 
attenuation features versus redshift. This will also enable to not only probe the 
present day EBL but to study its evolution with redshift. How such a study can be 
performed with the help of an NGCTS will be the topic of a second paper. 
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Appendix A. References for VHE spectral data 



Source name Redshift Reference 

MknSOl 0.034 Q7J 

PKS2155-304 0.116 [391. [361 

lES 1101-232 0.188 [51 

Table A.l: VHE spectra used in this publication. 
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